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Abstract. In this paper we present a scale free method to determine the cosmological parameters (fi m ,fiA)- The method 
is based on the requirement of isotropy of the distribution of orientations of cosmological filaments. The current structure 
formation paradigm predicts that the first structures to form are voids and filaments, causing a web-like structure of the matter 
distribution at high redshifts. Recent observational evidence suggests that the threads, or filaments, of the cosmic web most 
easily are mapped in Lya emission. We describe how such a 3D map can be used to constrain the cosmological parameters in 
a way which, contrary to most other cosmological tests, does not require the use of a standard rod or a standard candle. We 
perform detailed simulations in order to define the optimal survey parameters for the definition of an observing programme 
aimed to address this test, and to investigate how statistical and observational errors will influence the results. We conclude that 
observations should target filaments of comoving size 15 — 50 Mpc in the redshift range 2 — 4, and that each filament must be 
defined by at least four Lya emitters. Detection of 20 filaments will be sufficient to obtain a result, while 50 filaments will make 
it possible to place significant new constraints on the values of fi m and S1a permitted by the current supernova observations. In 
a future paper we study how robust these conclusions are to systematic velocities in the survey box. 
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1. Introduction 

Recent studies have successfully narrowed down the permit- 
ted parameter space for the cosmological parameters f2 m and 
£1a. Confidence intervals defined by observations of distant su- 
pernovae (Riess et al. 1998, Perlmutter et al. 1999) and those 
resulting from the high-resolution observations of the cosmic 
microwave background radiation (CMB) (Jaffe et al. 2000) 
meet almost orthogonally, and as a result they combine to 
bracket a domain of high probability. There remains, however, 
the question of possible systematic effects which may produce 
significant errors (e.g. Simonsen & Hannestad 1999; Rowan- 
Robinson 2002). It is therefore important that new and inde- 
pendent methods to determine f2 m and Q\ should be sought 
and exploited. 

The classic cosmological tests (Sandage 1961) involve the 
use of either standard rods or standard candles. The supernovae 
projects for example, are based on the use of standard candles, 
while recent proposals to use the distribution of Lyman Break 
Galaxies and arrangements of Lya-forest systems (Roukema 
2001 and McDonald & Miralda-Escude 1999 respectively) in- 
volve the use of standard rods. The potential worry with all 
such methods is that the hypothetical standard rods or candles 
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may in reality not be good standards because of evolutionary 
effects, or simply because the evolution of structure may be a 
fractal process without any preferred scale. 

The simplest and cleanest cosmological tests one might 
imagine would be a purely geometric test not involving mea- 
surements of standard rods. The first such method was pro- 
posed by Alcock & Paczyriski (1979), who considered an ide- 
alized set of objects distributed spherically symmetric and on 
average following the Hubble flow. The idea is beautiful in its 
simplicity, but the requirement of spherical symmetry of the 
distribution (requiring a length scale to be the same in all di- 
rections) is in reality an indirect use of a standard rod, even 
if at a single point in redshift space. The hypothetical spheri- 
cal distribution might as well be ellipsoidal, which would then 
invalidate the result. The question is if it is at all possible to 
devise a geometrical test not requiring any length scale at all. 

Numerical simulations of early structure formation based 
on the Cold Dark Matter scenario (Klypin & Shandarin 1983, 
White et al. 1987, Rauch et al. 1997) show that high redshift 
star-forming regions tend to align themselves in long filaments 
with a higher-than-average matter density. This has recently 
been confirmed observationally by M0ller & Fynbo (2001), 
who also pointed out that such structures offered three inde- 
pendent cosmological tests, all of which were purely geomet- 
ric, and none of which required the early structure formation to 
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have a preferred length scale. The only requirement is that of 
global isotropy. The first of those three tests concerns the dis- 
tribution of filament inclination angles which, in an isotropic 
universe, must be random. In this paper we present a detailed 
analysis, based on Monte Carlo simulations, of this first fila- 
ment test. 

The paper is organized as follows: In section 2 we introduce 
the method, in section 3 we describe the details of the simula- 
tions and interpret the results in a cosmological context to see 
which constraints we may put on the cosmological parameters, 
in section 4 we address some additional sources of statistical 
and observational errors, and use the results to define the op- 
timal survey parameters for an observing programme, and fi- 
nally in section 5 we summarize the results and conclusions. 
We have included a short appendix, which briefly introduces 
the notation and the necessary cosmological relations. 

2. F(n m , A , z): The line of sight scale factor 

The cosmological principle tells us that any class of elongated 
objects or structures will, at all redshifts, display a random dis- 
tribution of orientations of their major axes. Elongated struc- 
tures that are so large that they are everywhere anchored in 
the local Hubble flow are usually referred to as 'filaments' 
and 'walls', and for such objects the cosmological redshift will 
vary along the structures and thereby reveal their orientations 
in 3D redshift space. The basis of the cosmological test de- 
scribed here is to determine the sets of (f2 m , ^a) which, at any 
given redshift, will make the observed distribution of orienta- 
tions anchored in the Hubble flow conform to this requirement 
of isotropy. 

As a simple example let us first consider an area on the sky 
in the form of a square with sides of angular size 0. At any 
given redshift z this will correspond to a square in the plane of 
the sky with sides of proper size 



Wteue = T5-/w(^m, ^A, z) 
no 



(1) 



where prescriptions for calculation of /w (as well as of /l be- 
low) are provided in appendix A. If we further consider a red- 
shift interval Sz centred on z, we have defined a box of proper 
length 



itruc — ^A, z). 

no 



(2) 



Let the box be filled with a large number of filaments, and 
let their distribution of orientations be isotropic. If one were 
to either stretch or squeeze the box along the line of sight 
(along the redshift direction), the angular distribution of the fil- 
aments would (because they are anchored in the Hubble-flow) 
no longer remain isotropic. 

For illustration purposes it is useful to think of a box inside 
which an isotropic network of rubber bands has been mounted 
between opposing walls. If we stretch the box along one direc- 
tion the angles between the rubber bands and the stretching di- 
rection will all decrease. If, instead, we were to stretch the box 
by that same factor in all directions, the angles would remain 



unchanged. Hence, in considering observationally only the dis- 
tribution of angles, one is not gaining any information about 
true sizes. By the same token, since the true size of an object 
falls out of the equations, the test that we describe here does 
not require the existence of a 'standard rod' at any redshift. 

Because the actual length and width of the box is irrelevant 
for our test, the only numerical value of interest is the "change 
of scale" between the two directions (along the line of sight and 
perpendicular to the line of sight). What we are interested in is 
therefore the ratio between the length and the width 



L_ _ Sz f L (n m ,tt A ,z) 

w ~ cj> f w {n m ,n A ,z) 



(3) 



and how the apparent value of that ratio relates to its true value 



F{n m ,n A , z ) 



(W)ap 
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Note that Ho, <fi, and Sz have cancelled out of the ratio 
F(Q, m , Oaj z). In what follows we shall refer to this ratio as 
the "line of sight scale factor" . 

Returning briefly to the example above, for a given triple of 
observables (4>,z,Sz) and an assumed cosmology (Cl m ,Cl\), 
the box is fully defined. If one now assumes a different cosmol- 
ogy, the inferred size of the box will change. For most choices 
of cosmology, the change in length and width of the box will 
be different (hence the value of F will change), but it is always 
possible to find sets of (fi m , ^a) with the combined effect of 
changing both length and width by the same factor, thereby 
keeping F constant. In Fig. [j] we show a set of such "isoscale 
curves" (curves of constant F) for a redshift of 3. Our test can 
now be formulated simply: "Determine the sets of (f2 m , JIa) 
for which F = 1". The remainder of this paper is concerned 
with the formulation of how this is done in practice. In addition 
we use Monte Carlo simulated data to determine what level of 
accuracy may be reached with current telescopes and instru- 
mentation. 

3. Limits on F from Monte Carlo simulations 

Consider a large number of idealized, randomly aligned, lin- 
ear and thin filaments. Let the inclination angle 9 of a filament 
be the angle between the filament major axis and our line of 
sight to the filament. A filament pointing towards us will hence 
have inclination angle zero, and the cumulative distribution of 
isotropically distributed inclination angles can be shown to be 



P(x) = 1 - 



1 



VT 



(5) 



where x £ [0, tt/2] and P(x) is the probability of observing an 
inclination angle 9 < x. If we now use a wrong set of values 
(f2 m , f2 a), then the apparent distribution is no longer isotropic 
and the accumulated distribution of inclination angles is instead 
given by the general expression 



P{x,F) = 1 



1 



Vl + F 2 tan 2 x 
where F is the line of sight scale factor (Eq. 4). 



(6) 
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Contours of F at z = 3 



50 filaments with F=1 




Fig. 1. Solid: Contours of constant value of the dimensionless 
line of sight scale factor F (here for z — 3 and normalized to 
F = 1 for n m = 0.3, A = 0.7). Dotted: 68.3%, 95.4% and 
99.7% confidence regions obtained by supernovae observations 
(Riess et al. 1998). 



Inclination angle 



Fig. 2. Cumulative distribution of Monte Carlo simulation of 
50 filaments, compared to expected distributions for a range of 
values of the line of sight scale factor F. 



Precision on F for a number of filaments 



Given a set of observed inclination angles, a K-S test (Press 
et al. 1989) against P(x, F), can be used to determine the 
value of F which best fits the observed distribution. At present 
there is no data set available large enough to place useful con- 
straints on the cosmological parameters. Instead we have cre- 
ated Monte Carlo data in order to determine how large a data 
set we shall need to be able to set interesting limits. The K- 
S test is only valid for fairly large samples, so in addition we 
have tested a number of alternative estimators all based on the 
comparison of cumulative distributions. The classic K-S test 
only seeks to minimize the largest difference between the two 
cumulative distributions. For small samples we found that the 
scatter of the fit was smaller when we minimized the sum of 
the square of the difference at each entry in the sample (Eq. 7), 
while for large samples the two methods gave identical scatter. 

The samples of simulated filaments were computed by ran- 
domly selecting N inclination angles, 6i, such that the values of 
cos 0i formed a uniform distribution between and 1. For each 
such sample we used the extended K-S test described above 
and determined the best fitting value of F by minimization of 
D(F): 



N 



D(F)=J2 d * 



(7) 



i=i 



where di = max i 



P(0i,F)\ 



P(9 Z ,F)\), and P 



n * y" i > - 1 i\ ' I n 
is defined in Eq. 6. As an illustration we show in Fig. 2 the 

cumulative distribution of a realization of 50 filaments (his- 
togram) and the predicted distribution for a range of values of 
F (lines) based on Eq. 6. 

It is now interesting to ask how many filaments one would 
need to observe, to be able to place cosmologically interesting 
limits on the observed value of F. To address this question we 
repeated the simulations for the range 10 to 100 filaments in 
steps of 10, and in each case running 2 x 10 4 simulations. The 




Tiber of filaments 



Fig. 3. Plot of the uncertainty with which we can determine 
the stretching factor F as a function of how many filaments 
are in our sample. This is the result of 2 x 10 4 simulations for 
each number of filaments. The upper and lower 1, 2 and 3rr 
uncertainties are fitted with exponentials. The result is given in 
table [j] 



results are shown in Fig. ^| The one-sided 1, 2 and 3er curves 
were found as the 15.9%, 2.28% and 0.13% quantiles. The 1, 
2, and 3er curves are fitted well by the exponential functions 
given in Table [j]. The median is consistent with a straight line 
with intercept 1.00 and slope 0.000. 

As illustrated in Fig. 1, each value of F represents, at a 
given redshift, a curve in the (fi m ,fiA) diagram. Using the 
cosmological relations given in appendix we can therefore 
transform our limits on F (Fig. [}) into a set of curves confin- 
ing the permitted cosmologies (Fig. Q). Part of the parameter 
space in Fig. ^ represents "bouncing" universes that do not ex- 
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Table 1. The na errors fitted with F(N) = ae bN + c. 



na 


a 


b 


c 


+3cr 


2.67 


-0.0650 


1.45 


+ 2a 


1.25 


-0.0560 


1.27 


+1(7 


0.482 


-0.0507 


1.13 


-la 


-0.313 


-0.0451 


0.887 


-2a 


-0.498 


-0.0396 


0.791 


-3a 


-0.671 


-0.0417 


0.694 



pand monotonically from a Big Bang. This corresponds to the 
region in the figure where 51a > 4f2 m a; 3 , in which 

_ J cosh [| cosli^f^ 1 - 1)] when < £l m < \ 
X ~ { cos [| cos" 1 ^,- 1 - 1)] when § < Q m ' 

(Felten & Isaacman 1986). In this region of the figure, the ex- 
pression (1 + z) 2 (Cl m z + 1) — Qaz (z + 2) under the square- 
root in the equation for t| (see Appendix A) becomes negative 
for some z (the sign of the expression shifts when the universe 
changes from expansion to contraction and vice versa). We do 
not calculate F in this region. 

The contours of F in Fig. |] are calculated for z = 3, for 
20 and 50 filaments, and are normalized to F = 1 for an 
Sl m = 0.3, ft\ = 0.7 universe. For easy comparison, the con- 
fidence regions obtained by supernovae observations (Riess et 
al. 1998) are also plotted. Two interesting points are immedi- 
ately obvious from Fig. [I| First it is seen that the isoscale curves 
are mostly horizontal, and that they place tight upper limits on 
J7a almost independently of r2 m . This is contrary to the corre- 
sponding confidence limits from the supernovae observations 
which span diagonal regions in the diagram. Secondly it is seen 
that already with a sample of 20 filaments one would set upper 
limits on the value of tt\ intersecting those currently set by the 
high z SN studies. Such a limited study would therefore not 
only provide a new and independent way to determine 57 a, it 
would already serve to further constrain the parameter space 
allowed by the SN studies. For 50 filaments the intersections 
would be such that they would reduce the area of the 1 a upper 
confidence limit set by the SN studies by 30%. 

4. Observational implementation of the test 

4. 1. The optimal redshift range 

For the planning of an observing campaign to carry out the 
proposed test, it is useful first to consider at which redshifts 
the test is most sensitive to realistic values of 57a, and how one 
may best identify filaments at those redshifts. In Fig. || we plot 
F as a function of redshift for several different cosmological 
models, also here normalized so that F — 1 corresponds to the 
57 m = 0.3, 17a = 0.7 model. The dotted and dashed curves 
below the F = 1 line correspond to 57a = models with 
Q m = 1 and f7 m = 0.3 respectively. The dot-dashed curves 
above the F = 1 line correspond to models with 57a > 0.7 (see 
caption for details). All curves converge at F = 1 in the limit 
of zero redshift. It is seen that for all models with S7a < l,Fi& 
almost independent of redshift in the redshift range z « 1.5—4, 



Constraints on (H m ,n A ) for 20 filaments at z = 3 




1 2 

Constraints on (Q m ,U A ) for 50 filaments at z = 3 




12 3 

Fig. 4. Solid: Confidence limits on F for observations of 20 
(upper plot) and 50 (lower plot) filaments at z = 3, normalized 
to a„ = 0.3, ft A = 0.7. Dotted: The 68.3% (smallest), 95.4% 
and 99.7% (largest) confidence regions obtained by supernovae 
observations. 



while F grows monotonically with redshift for models with 
51a > 1- Hence, in order to test against any cosmology with 
57a < 1 one should select the redshift in the range 1.5 — 4 
where filaments are most easily detected. 

4.2. Filament identification and sparse sampling 
errors 

Lya narrow band imaging has proven an efficient technique 
for identification of high redshift filaments (M0ller and Warren 
1998, M0ller and Fynbo 2001). Confirming spectroscopy then 
maps out filaments as strings of separate star forming regions, 
each glowing in Lya. In what follows we shall refer to those 
Lya emitting regions using the shorter name 'LEGOs' (Lya 
Emitting Galaxy-building Objects). Since each filament in this 
way is mapped by a finite number of objects, the accuracy with 
which one may determine its inclination angle will be a func- 
tion of the number of LEGOs defining it. 

As for the optimal redshift range and number of filaments, 
it is useful to consider what may be the minimum number of 
LEGOs needed per filament for an adequate determination of 
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Fig. 5. We plot F as a function of redshift for a range of cosmo- 
logical models. F = 1 corresponds to the £l m = 0.3, fl\ = 0.7 
model. The dotted and dashed curves below the F = 1 line 
correspond to J1a = models with Q m = 1 and £l m = 0.3 re- 
spectively. The dot-dashed curves above the F = 1 line corre- 
spond to models with (O m ,Q A ) = (0.2,0.8), (0.1,0.9), (0.3,1.2), 
(0.3,1.4), (0.3,1.5). 

its inclination angle. To address this question we ran a se- 
ries of simulations in which we randomly placed N LEGOs 
inside filaments defined as cylinders of length L and diam- 
eter D. We then determined the "observed" orientation vec- 
tor of each filament, (defined as the best fitting straight line 
through the TV points), and the "observed" length-to-diameter 
ratio, ^obs/^obs-Q The observed set of orientation vectors de- 
fine a solid angle centred on the true orientation vector, from 
which the sparse sampling error (er sp ) may be determined di- 
rectly as the RMS of the 2D distribution of the observed set 
of intersections between the orientation vectors and the sur- 
face of the unit sphere. In Fig. 6 we show contours of con- 
stant cr sp as a function of N and the £ bs/-Dobs ratio. As one 
would expect er sp is small for long/thin filaments but larger for 
short/wide filaments. It is also seen that for filaments defined 
by fewer than four LEGOs cr sp grows rapidly. M0ller & Fynbo 
(2001) reported the detection of a filament at z — 3.04 defined 
by eight LEGOs. Assuming (O m ,f2A.) = (0-3, 0.7) and using 
the procedure outlined above, we find L Q b s = 20.2 Mpc and 
-Dobs = 3.2 Mpc (comoving) for the z = 3.04 filament, and 
an inclination angle of 25.5°. Its position is marked in Fig. 6 
and it is seen that the sparse sampling error on its orientation is 
re 1.9°. 

In order to propagate the sparse sampling errors shown in 
Fig. 6 to errors on the determination of the line of sight scale 
factor F, we repeated the Monte Carlo simulations of Sect. 3, 
now including random errors on filament orientations as de- 
scribed above. The goal was to determine how large intrinsic 

1 The length of the cylinder is found as L \, a = \/V2 RMS(z), and 
the diameter D ^ s = y8 RMS(r), where z is the position of a point 
projected onto the axis of symmetry, and r is the distance of a point to 
the axis of symmetry. 
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Fig. 6. Contours of constant sparse sampling error on filament 
orientation as a function of the number of LEGOs defining the 
filament and the length-to-diameter ratio L bs/-D bs- The con- 
tours are er sp = 0.5°, 1°, 2.5°, 5° and 10°, and the cross marks 
the filament detected by M0ller & Fynbo (2001). The contours 
are based on 10 6 runs at each point. 



errors on the filament orientations we could tolerate before the 
error analysis summarised in Fig. 3 would be seriously com- 
promised. The Monte Carlo simulations were repeated for 20, 
50, and 80 filaments, and for er sp in the range 0° to 10°. The re- 
sults are summarised in Table 2 where we list errors on F from 
Fig. 3, and from the full analysis including sparse sampling er- 
rors. It is seen that even for cr sp as large as 10° the additional 
errors are insignificant. From Fig. 6 we see that a sp < 10° al- 
ways is achieved for filaments with L bs/-D b s > 3 defined by 
four LEGOs or more. 

Table 2. Errors (Ict) on determination of the line of sight scale 
factor F (Fig. 3) including sparse sampling of individual fila- 
ments. The results are based on 2 x 10 5 runs. 







Filaments 






20 


50 


80 


0.0° 


0.259 


0.162 


0.128 


5.0° 


0.265 


0.166 


0.131 


10.0° 


0.283 


0.179 


0.140 



In conclusion of this section, a first rough determination of 
(f2 m , J7a) can be obtained with the discovery of 20 filaments 
each defined by four LEGOs, or 80 LEGOs in all. With 50 fil- 
aments, or 200 LEGOs total, the errors are such that they will 
place significant new constraints on the values of fl m and f2y\ 
permitted by the current supernova observations. 

4.3. Optimal survey box size 

Up till this point we have discussed the error budget caused 
by purely statistical errors only; errors arising from a finite fil- 
ament sample and from the finite number of objects defining 
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each filament. Those errors are scale invariant and place no 
constraints on the size of the survey box. For observational rea- 
sons there is an upper bound to the size of the volume we are 
able to survey to the necessary flux limit. On the other hand 
we can only identify filaments of sizes smaller than the cho- 
sen survey box, and if that box is too small then measurement 
errors and errors caused by motion relative to the Hubble flow 
could become dominant. One may consider two types of mo- 
tion relative to the Hubble flow; "random motion" and "bulk 
flow" (streaming). 

4.3.1 . Random velocities and measurement errors 

Random velocities and measurement errors are most efficiently 
treated together because errors on the redshift are indepen- 
dent from random radial motion and therefore adds in quadra- 
ture. Let us consider the filament reported by M0ller & Fynbo 
(2001). In this case the error on the redshift determination of 
the Lya lines was 50 km s _1 (Fynbo, M0ller & Thomsen 
2001). At z s=s the typical peculiar velocity is 200 km s _1 
(Branchini et al. 2001). This includes both random motion 
and systematic flows towards large dark matter haloes, and 
can therefore be considered a firm upper limit to the random 
peculiar velocities at high redshifts. Only the radial compo- 
nent (100 km s _1 ) of the random velocity vector is observed, 
and in combination with the typical redshift measurement er- 
ror we have a 112 km s _1 radial error. This corresponds to 
dz = 0.0015 at z = 3, or to 9% of the length of the filament 
(390 kpc in an £! m = 0.3, Qa = 0.7, h = 0.65 universe). To 
understand the effects on a survey for filaments, let us consider 
the two limiting cases of filaments with 9 — 0° and 9 = 90°. 
For 9 = 0° the filament is aligned along our line of sight, and 
the additional scatter (RMS = 9% of total length) along the fil- 
ament will simply make it look at bit longer than it is in reality. 
Hence the L b s /-D bs ratio will be overestimated by a minis- 
cule amount, and a sp will be similarly underestimated. Since 
a sp is already irrelevantly small, this will have no consequences 
for the cosmological test. Considering now the case 9 — 90° (a 
filament perpendicular to the line of sight) we see that in this 
case the measured diameter of the filament (800 kpc) will be 
increased (along the line of sight) by the added scatter, result- 
ing in an added error on the determination of the inclination 
angle. The added error can be found from Fig. 6 via calculation 
of the predicted £ bs/-Dobs ratio. The filament considered here 
would (had it been found to have 9 = 90°) have had a diame- 
ter increase of 70% along the line of sight. Hence L bs/Ajbs 
would have changed from 6.6 to 3.9, and the corresponding 
sparse sampling error from 1.9° to 3.4°. From Table 2 we see 
that this 70% increase has no effect, but with another factor of 
two those errors would become dominating. 

Errors perpendicular to the line of sight are typically a 
fraction of an arcsec, or about two orders of magnitude less 
than those along the line of sight. Those can be completely ig- 
nored. A side-effect of this "orders of magnitude difference" 
is that thin filaments with small inclination angles will ap- 
pear observationally as thin filaments (strings), while thin fila- 
ments with large inclination angles will appear to have cross- 



sections which are elongated along the line of sight (ribbons). 
Observational determination of the magnitude of this "ribbon 
effect" will make it possible to determine the actual RMS of 
the non-Hubble flow velocities at z = 3 (the 200 km s _1 we 
used in this example is, as mentioned above, to be regarded as 
an upper limit). 

The 20 Mpc filament considered above remains fully de- 
fined when we include expected random errors. However, we 
see from Fig. 6 that a similar filament defined by only 4 LEGOs 
would move close to the curve marking the 10° sparse sam- 
pling error. Therefore, because of the expected magnitude of 
the errors caused by random motion, it will in general not be 
possible to determine inclination angles with sufficient preci- 
sion for filaments significantly shorter than 15 Mpc. Filaments 
significantly larger than 50 Mpc are still prohibitively difficult 
to identify for observational reasons, so at present an optimal 
search for filaments at z = 3 must target filaments of comoving 
length 15 — 50 Mpc. The optimal survey box has sides of length 
50 - 60 Mpc. 

4.3.2. Systematic peculiar velocities 

The final point we shall consider is the effect of bulk flow in 
the survey box. First we address the scale dependence. Any 
bulk motion on the same scale as the filaments themselves, or 
larger, will simply result in a translation of an entire filament 
and hence have no influence on the distribution of angles. Bulk 
motion on scales much smaller than the filaments will appear 
as random motion (treated in the previous section). Therefore, 
only bulk motion in the regime between those two limits (sev- 
eral, but not all, LEGOs in a single filament) is relevant. 

It is at present not easy to predict precisely what sort of 
bulk motion will be present on this scale. In the local universe 
we see large scale streaming towards massive dark matter ha- 
los, but since this is the result of acceleration over a Hubble 
time one would expect similar streaming at higher redshifts to 
have smaller amplitude. In simulations of the early universe 
two types of flows are seen. First it is often seen that proto- 
galaxies in the filaments stream along the filaments towards the 
nodes, second it is seen that nodes are attracted to each other. 
The "streaming towards nodes" will have the effect of mak- 
ing the filaments appear longer than they are (stretching), the 
movement of nodes towards each other will have the opposite 
effect (squeezing). To first order one might expect that the two 
effects cancel each other, but it is equally possible that one of 
the two dominates. If the latter is the case, it will cause a sys- 
tematic error on the determination of F. 

In the previous sections we have outlined how to determine 
F from a set of data. Those data could be observational, but 
they could equally well be taken from a simulation. The only 
way to address the possible effect of intermediary scale bulk 
motion is to apply our method to sets of filaments in simulated 
volumes. We shall return to the results of this work in a separate 
paper. 
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5. Summary and conclusions 

A significant number of 8 — 10 meter class telescopes became 
operational during the past decade, thereby moving several cos- 
mological tests from the domain of speculation into the domain 
of observation. The currently most successful test relies heav- 
ily upon the assumption that supernovae of type la are standard 
candles. Other proposed tests, based on e.g. distances between 
Lyman Break galaxies, seek to exploit an expected standard 
rod. 

In this paper we have presented a detailed discussion of a 
test which relies on neither standard candles nor standard rods. 
The test makes use of only the requirement of isotropy, and 
of the prediction from numerical simulations of the high red- 
shift universe that it has a filamentary structure. Direct obser- 
vational evidence that such filaments may indeed be found was 
presented by M0ller & Fynbo (2001) who also proposed three 
cosmological tests based on the existence of such filaments. 
Here we have in detail investigated and described the imple- 
mentation of the first of these tests. We have used Monte Carlo 
simulations to determine the accuracy with which one may ob- 
tain values of ft m and and we have described how addi- 
tional statistical errors and observational effects will affect the 
results. Guided by the detailed Monte Carlo simulations, and 
taking into account the capabilities of current instrumentation, 
we have considered what would be the optimal survey parame- 
ters for an observing campaign aimed at an implementation of 
this test. 

We find that the constraints this test can set are mostly on 
ft\, or in other words the limiting curves in the usual (f2 m , 
diagram are mostly horizontal. This means that they in- 
tersect the probability curves from the z < 1 SN la projects. 
Already with a sample of only 20 filaments the SN la results 
will be intersected, but with a sample of 50 filaments the area 
of allowed values in the (il m , Qa) diagram will be significantly 
reduced. The optimal redshift range for the test is found to be 
z = 1.5 — 4. In the range z — 2 — 4 filaments can conve- 
niently be detected using deep searches for faint Lya emitters 
(LEGOs). Each filament to be used in this test must have an 
observed length-to-diameter ratio of no less than 3, and must 
be defined by no less than four objects. 

We have considered the effect of random and systematic 
non-Hubble flow velocities. There are no observational deter- 
minations of peculiar velocities at high redshifts yet, so their 
impact is discussed in qualitative terms mostly. A detailed anal- 
ysis of those effects, which are important as they set the mini- 
mum scale size of structures which will be useful for the test, 
will be presented in a future paper. However, using the local 
RMS value for peculiar velocities (which we consider a realis- 
tic upper limit), we find that the optimal survey box is a cube 
with sides of length 50 — 60 comoving Mpc. 

Three surveys are currently underway which should pro- 
duce data useful for this test: i) the Large Area Lyman Alpha 
(LALA) survey (Rhoads et al. 2000) that targets Lya emitters at 
redshift 4.5 and 5.7, ii) a survey conducted with the SUBARU 
telescope which has resulted in a large sample of candidate Lya 
emitters at z — 4.86 (Ouchi et al. 2002), and iii) in the oppo- 
site end of the optimal redshift window there is a survey on the 



Nordic Optical Telescope (NOT) to map out a slice of the z = 2 
universe in Lya emission (M0ller et al. in prep). The NOT sur- 
vey started in 2001 and takes advantage of the excellent UV ca- 
pabilities of the NOT instrumentation which was demonstrated 
in similar observations of fields containing z = 2 host galaxies 
of Gamma Ray Bursters (Fynbo et al. 2002). Currently there 
are no large area Lya surveys targeting z = 3, and there are 
no large area Lya surveys conducted on the ESO VLT. Smaller 
targeted programmes aimed at z — 3 fields have been success- 
fully conducted on the VLT (Fynbo et al. in prep), confirming 
that the necessary volume density of Lya emitters will indeed 
be found. 
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Appendix A: Cosmology with A / 

We follow the treatment given in Longair (1998), which uses 
the metric in the form: 

ds 2 = c 2 dt 2 ~R 2 (t)[dr 2 + TZ 2 sm(r/TZ) 2 (de 2 + sm{9) 2 d<j) 2 )],(AA) 

where 1Z is the radius of curvature of the universe, and R is 
the scale factor normalized to 1 at the present epoch. The radial 
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coordinate r is the metric distance at time t, and it can be found 
as a function of z by integrating the differential equation 

— = — {{l + z) 2 {Cl ln z + l)-n A z{z + 2)J , (A.2) 

where f2 m and f^A are given by the present mean density p 
and the cosmological constant A, respectively: 

SirGppc 2 Ac 2 

" m = ^F' " A = 3i?f (A3) 

The radius of curvature K can be determined from the relation 

1 O m + Q A - 1 

tv = vJh! • (A ' 4) 

Longair defines the distance measure 

Ksm(r/K) if l/K 2 > 
D = { r if l/K 2 = (A.5) 

Ksmh{r/K) if l/K 2 < 

from which the angular diameter distance Da can be derived 
as Da = D/(l + z). The proper radial distance is determined 
from the differential (£) = + z). 

\ az / prop az ' ^ ' 

We can now finally derive the two functions and /w 
introduced in Sect. 2. As Wtruc = 4>Da the function is 
simply given by H Dj\. Similarly, the proper length L true of 
a box defined by two redshifts z min and z max (5z = z max — 

^min < (z max + z min )/2)isgivenby (fjj) prop hence /l = 

#o (tt) 

v \dz / prop 



